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ABSTRACT 

We present the analysis of a sub-Damped Lyman-o- system with neutral hydrogen column density, log N(H°) (cm -2 ) = 20.0+0.15 
at z abs = 2.69 toward SDSS J1237 14.60+064759.5 (z cm = 2.78). Using the VLT/UVES and X-shooter spectrographs, we detect H 2 , 
HD and CO molecules in absorption with logJV(H 2 , HD, CO) (cm 2 ) = 19.21+° J|, 14.48±0.05 and 14.17±0.09 respectively. The 
overall metallicity of the system is super-solar ([Zn/H] = +0.34 relative to solar) and iron is highly depleted ([Fe/Zn] = -1.39), 
revealing metal-rich and dusty gas. Three H 2 velocity cxomponents spanning ~ 125 km s" 1 are detected. The strongest H 2 component, 
at Zabs = 2.68955, with log N(H 2 ) = 19.20, does not coincide with the centre of the Hi absorption. This implies that the molecular 
fraction in this component, fm = 2A'(H 2 )/(2A'(H 2 )+A'(H )), is larger than the mean molecular fraction (fm) = 1/4 in the system. This 
is supported by the detection of Cl° associated with this H 2 component having N(C\°)/N(C\ + ) > 0.4. Since Cl° is tied up to H 2 by 
charge exchange reactions, this means that the molecular fraction in this component is not far from unity. The kinetic temperature 
derived from the J = and 1 rotational levels of H 2 is T = 108+j* K and the particle density derived from the C° ground-state fine 
structure level populations is «ho ~ 50-60 cm 4 . We derive an electronic density <2 cm" 3 for a UV field similar to the Galactic one and 
show that the carbon to sulphur ratio in the cloud is close to the solar ratio. The size of the molecular cloud is probably smaller than 
1 pc. Both the CO/H 2 = 10~ 5 and CO/C° ~ 1 ratios for fm > 0.24 indicate that the cloud classifies as translucent, i.e., a regime where 
carbon is found both in atomic and molecular form. The corresponding extinction, A v = 0.14, albeit lower than the definition of a 
translucent sightline (based on extinction properties), is high for the observed H° column density. This means that intervening clouds 
with similar local properties but with larger column densities (i.e. larger physical extent) could be missed by current magnitude-limited 
QSO surveys. The excitation of CO is dominated by radiative interaction with the Cosmic Microwave Background Radiation (CMBR) 
and we derive T ex (CO) = 10.5+^ K when 7" C mbr(z=2.69) = 10.05 K is expected. We measure W(HD)/2A'(H 2 ) = 10" 5 . This is about 
10 times higher than what is measured in the Galactic ISM for fm = 1/4 but similar to what is measured in the Galactic ISM for 
larger molecular fractions. The astration factor of deuterium - with respect to the primordial D/H ratio - is only about 3. This can be 
the consequence of accretion of unprocessed gas from the intergalactic medium onto the associated galaxy. In the future, it will be 
possible to search efficiently for molecular-rich DLAs/sub-DLAs with X-shooter but detailed studies of the physical state of the gas 
will still need UVES observations. 

Key words. Cosmology: Observations - Galaxies: ISM - Quasars: Absorption lines - Quasars: Individual: 
SDSSJ123714.60+064759.5 



1. Introduction 

Studies of the Interstellar Medium (ISM) in the local Universe 
have shown that the neutral ISM presents a complex struc- 
ture, with cold and dense clouds immersed in a warmer and 
more diffuse medium. These different ISM phases should be 

Send offprint requests to: P. Noterdaeme 

* Based on observations carried out with X-shooter and the 
Ultraviolet and Visual Echelle Spectrograph (UVES), both mounted 
on the European Southern Observatory Very Large Telescope Unit 2 
- Kueyen, under Program IDs 082.A-0544(A), 083A-0454(A) and 
084.A-0699(A). 



detectable at high redshift by their absorption signatures in 
Damped Lyman-a ( DLA) systems observed in quasar spectra 
(iPetitiean et al.L[T992h . However, although t here are evidences o f 
the multiphase nature of DLA systems (e.g. lWolfe et all [2004). 
most of the intervening DLAs probe only warm (T ^ 3000 K) 
and diffuse (»h < 1 cm" 2 ) atom ic gas (e.g. IPetitiean et all 
2000: lKanekar & ChengaluiU 2003>. The reason is that the cross- 
sections of the different phases are quite different and it is not 
possible to sample them equally well. 

Searching for molecular hydrog e n in high redshift DLAs 
dLedoux et all 120031: IPetitiean et all l2006t iNoterdaeme et all 
2008a) is an efficient way of detecting colder and denser 
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neutral gas and to probe its physical conditions (e.g . 



Reim ers et al.ll2003HCui et all 120051: [Hirashita & FerraraLr2 005 ; 
Srianand et all 120051; iLedoux et all l2006bt iNoterdaeme et all 
2007a, b). These studies have shown that molecular hydrogen 
is confined in small clouds (pc-sized) with densities n ~ 1- 
100 cirT 3 and temperatures T ~ 70-200 K. The filling fac- 
tor of Fb-bearing clouds in DLAs is much less than one and 
only 10 % of the lines of sight through a DLA galaxy do inter- 
cept H2-bearing cl o uds do wn to a limit of N(H.2) ~ 10 14 cm" 2 
dNoterdaeme et all l2008ah . H2-bearing clouds in DLAs have 
small physical extents. Direct evidence for this is given by 
the fact that the intervening H2-bearing gas does not com- 
pletely fill th e beam from the broad line region o f the quasar 
Q 1232+082 dlvanchik et all 120101: iBalashev et all submitted). 
Nonetheless, the molecular fraction in DLAs remains small and 
typical of what is seen in Galactic diffuse atomic gas with 
fm = 2 N(H 2 )/(2Ar(H 2 ) + MH°)) < 0.1 and often muc h lower 
tha n this dLedoux et all 120031 INoterdaeme et all l2008al) . 

ISnow & McCalll d2006l) have classified Galactic interstellar 
clouds into the following categories: (i) diffuse atomic, with low 
molecular fractions; (ii) diffuse molecular, where the fraction of 
hydrogen in molecules becomes substantial (/h 2 > 0.1) but car- 
bon is still ma inly in ionised form (C + ); (i ii) translucent (first 
introduced by Ivan Dishoeck & Blackl 119891) . where the carbon 
makes the transition to molecular; and (iv) dense molecular, 
where both hydrogen and carbon are fully molecular. As dis- 
cussed above, most of the H2-bearing DLAs detected so far are 
part of the first, and maybe for some of them, part of the second 
categories. The fourth category may be difficult to detect in ab- 
sorption because of the high extinction such a cloud produces on 
the background source. 

Despite their highly interesting chemistry and their close 
connection with star formation, we know very little about 
translucent clouds (i.e., the third category) at high redshift. The 
small cross-section of these clouds and/or the induced extinc- 
tion of the light from the background sources can probably ex- 
plain the absence of detection in more than three decades of QSO 
absorption-line research. However, observing molecular-rich gas 
in absorption should be possible by selecting sightlines passing 
through or starting from star-forming regions. 

Since long-duration Gamma-Ray Bursts (GRBs) are known 
to occur within star-forming regions, absorption lines at the 
host-galaxy redshifts which are impri nted in GRB optical af- 
terglow spectra (e.g. jFvnbo et all 120091) are obvious targets to- 
wards this goal. Nevertheless, current samp les are characterised 
by a general lack of H2 detection (e.g., iFvnbo et all [2006; 
Tumlinson et afl. l2007l) . This is probably due to the still limited 
sample sizes as well as a bias against dusty - molecu lar-rich - 
lines of sight dLedoux et all 120091: IFvnbo et all 120091) . The first 
detection of both H2 and CO in th e low-resolution spectru m of a 
highly reddened GRB afterglow (Prochaska et all l2Q09h seems 
to confirm this scenario. Moreover, the observed molecular ex- 
citation is high in this case, indicating strong UV pumping from 
the GRB afterglow itself. 

With the large number of quasar spectra available in the 
Sloan Digital Sky Survey (SDSS), it becomes possible to se- 
lect the rare sightlines passing through intervening molecular- 
rich gas. However, due to the small cross-section of such clouds, 
an efficient selection must be applied. In the local ISM, car- 
bon is found to transition from a ionised state (C + ) to neu- 
tral (C°) and molecular form (CO) fr om the most superficial 
to the deepest par ts of the clouds (e.g. ISnow & McCalll [2006; 
BurghetafllSoiO). Fro m our Very Large Telescope (VLT) sur- 
vey for H 2 in DLAs dLedoux et all |2003; N oterdaeme et all 



2008a), it appears that C° is generally observed in the same com- 
ponents as H2. This is due to the photo-ionisation potential of 
C° being similar to the energy of photons that dissociate H2. 
However, the neutral fraction of carbon is generally small, prob- 
ably because the gas is not completely shielded. This explains 
the non-detection of CO in th ese H2-bearing DLAs , even down 
to iV(CO) ~ 10 12 cm- 2 (e.g. iPetitiean et all [2002h . Searching 
for systems with large column densities of neutral carbon could 
be an efficient way to select more shielded gas where other 
molecules can survive, without relying on a pre-selection based 
on the H° column density (i.e. the absorbers need not be DLAs). 
Since several C 1 lines are located redwards of the Lyman-o' for- 
est, it is possible to search for strong C 1 absorptions directly in 
SDSS spectra using automatic procedures. We therefore initi- 
ated a program to survey with the VLT such specific sightlines. 
Our selection has been very successful and already allowed us 
to detect carbon mono xide along QSO sig htlines for the first 
dSrianand et all 1200 8b) and second times ( Noterd aeme et all 
2009a]). We present here the third detection of CO, at z = 2.69 to- 
wards SDSSJ123714.60+064759.5(Zem = 2.78, hereafter called 
J 1237+0647). This is a beautiful and peculiar case for which de- 
tailed analysis of the physical properties of the gas is possible. 
We present the observations in Sect. [2] the measurement and re- 
sults in Sect.[3]and provide some discussion in Sect. [4] We sum- 
marise our findings in Sect. [5] 

2. Observations 

2.1. X-shooter observations 

We are conducting an observing campaign with X-shooter 
mounted on the Cassegrain focus of the VLT Unit 2-Kueyen tele- 
scope to study the molecular content of our complete sample of 
C° absorbers. As a test case for the sensitivity of X-shooter in 
the blue, we observed J 1237+0647 (g = 19.2) twice in service 
mode on February 24 (airmass 1.2; seeing 1.4") and March 3, 
2010 (airmass 1.3; seeing 1.2"), using a slit width of 1" in the 
UVB arm. Each observation run consisted in 1 h exposure taken 
in staring mode (see Table [TJ. This yields the nominal resolu- 
tion power of R = 5100 in the UVB arm and a signal-to-noise 
ratio of about 50 at ~500 nm. Data were reduce d using version 
0.9.5 of the preliminary ESO X-shooter pipeline dGoldoni et all 
2006) and the appropriate calibration data. The two individual 
spectra were then combined weighting each pixel by the inverse 
of the error variance. A portion of the X-shooter spectrum fea- 
turing CO is shown on Fig. [1] where several electronic bands 
of CO are clearly detected. These bands are resolved into in- 
dividual rotational levels in the UVES spectrum obtained with 
8.5 h of exposure time (inset figures). From this, it is apparent 
that X-shooter is the most efficient instrument to survey a com- 
plete sample of candidates down to quasar magnitudes as faint as 
r ~ 21.5, whereas using UVES would be excessively time con- 
suming. Then, but only in the case of detection, higher spectral 
resolution is needed to make a detailed analysis of the physical 
state of the gas as done in the following. 

2.2. UVES observations 

The quasar J 1237+0647 was observed in visitor mode on 
March 27 and 29, 2009 and April 27, 200 9 with the Ultraviole t 
and Visual Echelle Spectrograph (UVES: iDekker et all 12000). 
mounted at the Nasmyth B focus of VLT-UT2. The total expo- 
sure time on source is 8.5 h (see TableQ]). We used two dichroic 
settings (4x5400 s with 390+564 and 2x4500 s with 390+775) 
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Table 1. Journal of observations 



Instrument 


Date 


Setting 


Exposure time 


Resolving power 


SNR* 


UVES 
UVES 
UVES 


27-03-2009 
29-03-2009 
27-04-2009 


390+564 
390+564 
390+775 


2x5400 s > 
2x5400 s 
2x4500 s 


I 390B:51400 - 564L:50800 - 564U:49500 - 
" 775L:50200 -775U:48600 


10-40 


X-shooter 


24-02-2010 


UVB 


3600 s } 


5100 


50 


X-shooter 


03-03-2010 


UVB 


3600 s 



Notes. (a) 'B', 'L' and 'U' stand for respectively blue, lower red and upper red CCD. ( ' per pixel 




5250 5300 5350 5400 5450 

Observed wavelength (A) 

Fig. 1. Portion of X-shooter spectrum in the UVB. The inset fig- 
ures show 5 A-wide portions of the UVES spectrum around the 
position of the detected electronic bands of CO. 

to cover the wavelength range 3300-9600 A with small gaps at 
4517-4621, 5597-5677 and 7764-7809 A. The CCD pixels were 
binned 2x2 and the slit width adjusted to 1", yielding a re- 
solving power of ~50000 under seeing conditions of 0.9-1". 
Individual science spectra were reduced using the ESO UVES 
pipeline, which performs accurate sky subtraction while remov- 
ing cosmic ray impacts at the same time. The spectra were then 
combined using a dedicated IDL routine by weighting each pixel 
by the inverse of the error variance in that pixel and clipping 
residual cosmic rays impacts that remained after the cleaning of 
2D spectra. 

3. Analysis 

The system at z = 2.69 towards J 1237+0647 features numer- 
ous absorption lines from atomic (H°, O , C°, Mg°, Cl° and 
S°), singly-ionised (Fe + , Si + , Zn + , Ni + , S + , C + ), and molecular 
species (two isotopomers of molecular hydrogen: H2 and HD; as 
well as carbon monoxide: CO). 

We analysed the UVES spectrum using standard Voigt pro- 
file fitting techniques. The fits were per formed through y 2 - 
minim isation using the code FITLYMAN (Fontana & Ballester, 
1 1995b which is available as a context of the ESO-MIDAS 
data analysis software. The spectrum was normalised in the 
wavelength ranges of interest by fitting spline functions to 
regio ns free f rom absorption lines. Atomic data were taken 
from iMortonl d2003l) for metal lines, unless otherwise spec- 
ified. Wavelength s and oscill ator strengths were t aken from 



iMorton & Noreaul d!994l) and jp idelsberg & Rostasl (120031) for 
CO and from lAbgrall & Roueffl (120061) for HP. Updated wave 



lengths of H? Ly man and Werner bands were taken from 
Baill v et al] (|2010), with oscillator streng ths from the Meudo n 
grourn, based on calculations described in Abgr alletal.l(fl99l . 



http : //amrel . obspm.fr/molat/ 



Photos pheric solar abundances are taken from lAsplund et alJ 
J2009I) . 



3.1. Atomic hydrogen 

From the damped Lyman-o- absorption line (see Fig. |2}, we 
measure the total column density of atomic hydrogen to be 
log Af(H°)(crrr 2 ) = 20.00 + 0.15, whi ch is in agreement with 
the value measured automatically bv lNoterdaeme et all ((2009b) 
from the low resolution SDSS spectrum (20.15 + 0.28). The 
centroid of the Hi profile is well constrained by the Lyman- 
/3 and Lyman-y absorption lines. The large Doppler parameter 
(b ~ 100 kms _I ) required to fit the Lyman-/? and Lyman-y 
lines is a consequence of the presence of multiple components 
as testified by the clumpy profile of the O 1/11302 absorption 
line spread over ~350 km s _1 (see Fig. We recall that O 
closely follows H° because of favourable charge-exchange re- 
action. Unfortunately, because of strong saturation and blending 
effects, it is not possible to derive column densities in individual 
components and only the total H° column density along the line 
of sight is accessible. 

In Fig. [2] and subsequent figures and tables, the zero of the 
velocity scale is taken at the position of the CO component 
fobs = 2.68957, see Sect. 13.5b and the centroid of the three de- 
tected H 2 components at z abs = 2.68801, 2.68868 and 2.68955 
(Av = -127, -73, -1.6 km s _1 , see Sect. 13.31 1 are indicated by 
short vertical marks. Interestingly, the centroid of the atomic 
hydrogen absorption profile (vertical dotted line in Fig. [2] at 
z a b s = 2.69063) is shifted by about +86+10 km s _1 relative to 
the CO absorption feature. This, the clumpy Oi profile, and the 
large value of the ^-parameter of H 1 lines, all indicate that a 
significant fraction of the atomic gas is not associated with the 
molecular gas. We will discuss this further down in more details. 



3.2. Metal content 

3.2.1. Singly ionised species 

Absorption lines from detected low ionisation species are spread 
over about 350 km s _1 around the strongest component, which is 
also the component where CO and HD are detected. We used 
non-saturated transitions to derive accurate column densities for 
Fe + , Ni + , S + and Zn + . Lines from other species (C + , O , N°) are 
heavily saturated, preventing us to derive any meaningful value 
of the corresponding column densities. It is however possible 
to perform an accurate measurement of the Si + column density 
from the simultaneous use of Si 11 A 1808, which is very weak 
(below the 3 cr detection level), and Siii/11304 which is close to 
saturation. Note that although the strong Siii/11304 line reveals 
the presence of additional weak components, their contribution 
to the overall column density is negligible. The result of the 
Voigt-profile fitting is shown on Fig. [3] with the corresponding 
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Fig. 2. Measurement of the total column density of neutral 
atomic hydrogen at z — 2.69 towards J 1237+0647. The solid 
line represents the best one-component fit to the data. The 
centroid of the profile is indicated by a dotted vertical line. 
Uncertainties on the column density measurement are illustrated 
by the shaded area. The short tick marks indicate the position of 
the three H2-bearing components, the reddest of which also fea- 
tures CO and HD absorptions. The origin of the velocity scale, 
for this figure and all following ones, is defined at the position 
of the CO-bearing component at z a bs = 2.68957. 



parameters in Table|2] Both the S H/11253 and S n/11259 profiles 
are partially blended. Unfortunately, S n/11250 falls in a gap of 
the spectrum. Therefore, we provide only upper-limits for the 
blended components in Table [2] 3-cr upper-limits on the column 
densities of Ni + and Zn + for the undetected components are pro- 
vided in the table. Finally, we measure log N(Cr + ) < 12.2 at 3 cr 
from the non-detection of Cr H/12056. 

Total column densities and corresponding mean metallicities 
in the gas relative to solar are given in Table [3] Note that we 
do not apply ionisation correction since the presence of neu- 
tral and molecular species in the strongest components indi- 
cates the effect of ionisation on the overall abundances should be 
negligible. Indeed, even in the general population of absorbers, 
the ionisati on correction is on ly about 0.1 dex for A^(H°) = 
10 20 cm -2 dPeroux et al.L 120071) . The metallicity is super-solar 
with [Zn/H] = +0.34 and [S/H] = +0.15. Other species are 
depleted ([Fe/Zn] = -1.39, [Si/Zn] = -0.82). This indicates 
that a significant fraction of the refractory species is locked in 
solid phase onto dust grains. The presence of dust and its conse- 
quences are discussed in more details in Sect. 14.21 

3.2.2. Neutral carbon 

As the ionisation potential of neutral carbon (C°) is similar to the 
energy of the photons t hat destroy H?, C i is u sually a good tracer 
of the presence of H2 dSrianand et all [20051) . The expected po- 
sitions of several C 1 lines usually fall out of the Lyman-a for- 



Table 3. Summary of overall gas-phase abundances 



Species log A' (cm 2 ) mean abundance" 



H° 


20.00±0.15 




H 2 


19 21 +013 


</ H2 > = 0.24 


CO 


14.20+0.09 


log CO/H = -5.92 


Zn + 


13.02±0.02 


[Zn/H] = +0.34+0.12* 


S + 


15.39±0.06 


[S/H] = +0.15+0.13* 


Fe + 


14.57±0.01 


[Fe/H] =-1.05+0.12* 


Si + 


15.15±0.02 


[Si/H] =-0.48+0.12* 


Ni + 


13.48+0.03 


[Ni/H] =-0.86+0.12* 



Notes. {a) Abundances are given considering the total neutral hydrogen 
column density N(H) =N(Yl Q ) + 2N(H 2 ). ( *' Relative to solar abundances 
(Asplu nd et alll2009l) . 



est. We therefore initiated a program to search for molecules 
along QSOs selected upon the presence of Ci, as seen in the 
low resolution SDSS spectra. Because of this selection, it is not 
surprising to detect strong Ci lines in the UVES spectrum of 
J 1237+0647. The profile of C 1 absorption lines is complex and 
results from the blending of absorption lines from different com- 
ponents seen in different excitation levels (ground state: 3 Po, first 
excited level: 3 fi and second excited level: ^Pi). Nevertheless, 
the high signal-to-noise and high spectral resolution allow us to 
clearly identify eight components. Most of them are also de- 
tected in the first excited level, while only the strongest two 
are detected in the second excited level. The fit to C 1 lines is 
shown on Fig. [4] with the corresponding parameters given in 
Table[4] We considered all optically thin absorption lines but did 
not include in the fit weak absorption lines in the region around 
C 1/11277 where the placement of the continuum was uncertain. 
The main uncertainty in determining the C° column densities 
comes from the uncertainties on th e oscillator strengths. As in 
previous works from o ur group (e.g. Noterdaeme et al.Ll2007allbt 
Srianand et al ., 2008b) and ot hers in the fiel d dJorgenson et all 
2009), we used /- values fromlMortonl(l2003l) . Using the oscilla- 
tor strengths from lJenkins & Trippl d200lh results in 2 to 3 times 
lower column densities. 

The relative populations of the fine-structure levels of neu- 
tral carbon depend on the gas pressure. Since the kinetic tem- 
perature of the gas can be derived from the the relative pop- 
ulations of the low rotational levels of H2 (see Sect. 13.3b . 
it is possible t o measure the volum ic densi ty of the gas . 
From figure 2 of ISilva & Viegasll2002l (see also ISrianand et all 
2000) and taking into account excitation by collisions and 
by the Cosmic Microwave Background radiation, we can see 
that the measured ratios log N(C ,J=1)/N(C°,J=0) = -0.2 and 
log Af(C°,J=2)/Af(C°,J=0) = -1 in the main component, coincid- 
ing with the position of the CO component, correspond to a volu- 
mic density of the order of n H o ~ 50-60 crrT 3 for T ~ 1 10 K (see 
Sect. 13.31 ). Other components have similar fine-structure ratios, 
which indicate similar thermal pressure. The kinetic temperature 
is probably larger in all other components (which is verified at 
least for the two other H2-bearing components, see Tabled, im- 
plying smaller densities (e.g. n H ° ~ 1-10 ctrT 3 in the component 
at v = -127 km s" 1 ). 

3.2.3. Neutral sulphur 

The first ionisation potential of sulphur being of 10.36 eV, 
this element is hence usually observed in its first ionised 
state (S + ) in Damped Lyman-ff systems. In turn, sulphur is 
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Fig. 3. Fit to metal lines. The origin of the velocity scale is set at the redshift of the CO absorption (z = 2.68957). The positions of 
the three H2 components are indicated by short tick marks. The absorption seen at v = +55 km s in the Zn H/12026 panel (in blue) 
is due to Mg 1/12026 while that at v = -110 kms 1 on the Sin/11808 panel is due to SiA1807. Note that both S 11 lines are affected 
by blends. Finally Si H/11304 reveals the presence of additional weak components that are not fitted (see text). 



expected to be found in neutral form (S°) inside molecular 
clouds, where the surrounding UV field has been strongly at- 
tenuated. To date, S° has been detected only in QSO ab- 
sorbers where CO absorption is seen as well: in the sys- 
tems at z=2.4 2 towards SPSS J14391 2.04+1 11740.5 (hereafter 
J 1439+1117; ISrianand et all 120 08b) and at z =1.64 towards 
SDSS J160457.50+220300.5 (iNoterdaeme et all 120091 . 

Absorption lines from eight transitions of neutral sulphur are 
detected in the UVES spectrum of J 1237+0647. This strength- 
ens the claim that the presence of neutral sulphur flags molec- 
ular gas. Just like th e presence of neutral carbon is a good in- 
dicator of that of H2 dSrianan d et al.ll2005[. albeit with generally 
modest molecular fractions, INoterdaeme et al.ll2008al) . S 1 lines 
might well indicate the presence of CO. This is however of little 
practical use for pre-selecting CO-bearing DLAs from the low 
resolution SDSS spectra, since S 1 and CO lines are located in 
the same spectral region and have similar strengths. 

We used all detected S 1 lines to constrain the column density 
and b parameter. Two components are needed to properly fit the 
data, with resulting xl - 1 ■ The fr-value obtained is less than 
1 kms 1 for the main component, i.e., well below the spectral 
resolution (~6 km s ). However, b is well constrained, thanks to 
the relatively large range spanned by the oscillator strength val- 



Table 5. S° column densities 



2abs 


Av" (kms- 1 ) 


logA^S ) (cm- 2 ) 


b (km s 1 ) 


2.68953 


-3 


12.29+0.06 


1.1+0.5 


2.68958 


+ 1 


13.19+0.04 


0.7+0.1 



Notes. (n) Relative to z abs = 2.68957. 



ues. However, the fit is sensitive to the exact value of the spectral 
resolution. Therefore, to add confidence to the b and N measure- 
ments, we built the curve of growth for the detected S 1 lines, 
which does not depend on the spectral resolution (see Fig. |6j. 
The error on the equivalent width measurements are conserva- 
tive and take into account uncertainties in the continuum place- 
ment. From this figure, we confirm the small Doppler parameter. 
The measured column density nicely matches the sum of the in- 
dividual column densities in the two components derived from 
the Voigt-profile fitting. 
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Table 2. Column densities of metal species 



Zabs 


Av"(kms-') 


b (km s _1 ) 


Fe + 


Ni + 


\ogN>> 
S + 


(cm" 2 ) 
Zn + 


Si + 


Mg° 


2.68803 


-125 


5.2+0.2 


13.03+0.06 


12.18+0.13 


14.03+0.02 


11.46+0.15 


13.82+0.02 


<12.70 


2.68870 


-71 


4.6+0.1 


13.56+0.04 


12.62+0.05 


14.04+0.02 


11.86+0.06 


14.37+0.04 


<12.70 


2.68886 


-57 


4.8+0.3 


13.09+0.05 


<12.15 


13.83+0.02 


11.50+0.13 


<13.91 


<12.70 


2.68909 


-39 


6.2+0.5 


13.30+0.04 


12.07+0.19 


14.18+0.01 


11.86+0.06 


< 14.08 


< 12.70 


2.68925 


-26 


6.0+0.3 


13.37+0.04 


12.52+0.07 


14.01+0.02 


11.76+0.08 


14.04+0.09 


< 12.70 


2.68958 


+ 1 


14.5+1.5 


13.41+0.05 


<12.15 


14.17+0.07 


<11.55 


13.94+0.04 


<12.85 


2.68961 


+3 


5.4+0.2 


13.74+0.02 


12.79+0.04 


14.91+0.02 


12.75+0.02 


14.04+0.09 


13.03+0.02 


2.69024 


+55 


3.5+0.3 


12.90+0.06 


<12.15 


<13.50 


< 1 1 .40 


13.21+0.08 


<12.70 


2.69044 


+71 


14.8+2.0 


13.69+0.02 


12.67+0.09 


< 14.70 


<11.55 


14.37+0.01 


<12.85 


2.69086 


+ 105 


9.3+2.0 


13.46+0.03 


<12.20 


<14.80 


11.40+0.32 


14.09+0.02 


<12.75 


2.69099 


+ 115 


3.7+0.2 


13.75+0.02 


12.65+0.04 


<14.70 


12.12+0.04 


14.08+0.14 


<12.70 


2.69124 


+136 


3.4+0.5 


12.81+0.04 


<12.15 


<13.80 


11.15+0.29 


13.58+0.04 


<12.70 


2.69189 


+ 188 


19.5+1.7 


13.40+0.03 


<12.40 


<13.40 


<11.60 


14.04+0.05 


<12.95 


2.69220 


+214 


6.5+0.5 


13.12+0.03 


<12.15 


13.78+0.04 


< 1 1 .40 


13.78+0.02 


<12.70 



Notes. <n) In all figures and tables, the velocity is given with respect to the redshift of the CO component at z = 2.68957. ib) Upper-limits due to 
blends or saturated lines are indicated by '<', while 3 a upper-limits from non-detections are indicated by '<'. 



Table 4. Column densities of neutral carbon in fine-structure levels 





Av" 
(kms- 1 ) 


b 

(km s" 1 ) 


3 Po 


logJV(C°) (cm" 2 ) 


total 


2.68802 


-126 


1.2+0.1 


13.24+0.04 


12.80+0.03 


13.37+0.03 


2.68872 


-69 


11.8+0.6 


12.84+0.03 


12.70+0.05 


13.08+0.02 


2.68888 


-56 


3.3+1.3 


12.37+0.08 




12.37+0.08 


2.68906 


-41 


3.0+0.4 


12.82+0.02 


12.56+0.04 


13.01+0.02 


2.68926 


-25 


7.4+0.4 


12.94+0.02 


12.87+0.02 


13.21+0.01 


2.68956 


-1 


1.4+0.1 


14.67+0.04 


14.46+0.03 13.64+0.02 


14.90+0.03 


2.68960 


+2 


8.8+0.1 


13.98+0.01 


13.79+0.01 12.98+0.03 


14.22+0.01 


2.69097 


+ 114 


3.6+0.5 


12.70+0.02 


12.42+0.08 


12.88+0.03 



Notes. (fl) Relative to z abs = 2.68957. 



3.2.4. Neutral chlorine 

Chlorine, with an ionisation potential of 12.97 eV is a unique 
species among those that can be photoionised by photons with 
energy hv < 13.6 eV. The dominant form of chlorine is Cl + when 
hydrogen is mostly in the atomic form. However, neutral chlo- 
rine (Cl°) results from rapid exothermic ion-molecule reaction 
between si ngly-ionised chlorine (Cl + ) and H2 when H2 is opti- 
cally thick (ESGlZl). Therefore, the presence of neutral chlo- 
rine in the ISM is expected to be a good indicator of the presence 
of molecular gas. 

CI is clearly detected at z a t, s = 2.689560, associated with the 
strongest H2 component. We measure log/V(Cl°) = 13.01+0.02 
from the fit to the CI 1/11347 absorption line, with b - 4.5 + 
0.4 kms 1 (see Fig [7j. From the non-detection of Clii/11071, 
we derive log A^(C1 + ) < 13.4 at the 3 <x confidence level, which 
translates to f a o = N(Cl°)/(N(Cl°) + N(Cl + )) > 0.3. As the 
fraction of chlorine in neutral form is expected t o follow ap- 
proxi mately t hat of hydrogen in molec ular form (lJura & York 
1978, see also Sonnentrucker et al. 2002), the lower limit on / cl o 
indicates that hydrogen could be mostly molecular at the place 
where we detect Cl°. We indeed show in the next Section that 
the molecular fraction is particularly high in this component. 

3.3. Molecular hydrogen 

Molecular hydrogen is detected in three components at 
z abs =2.68801, 2.68868 and 2.68955, spread over - 12510ns- 1 . 



The strongest component at z a bs = 2.68955 also features HD 
and CO absorption lines. The UVES spectrum covers numerous 
Lyman bands (B'2,/ <— X'S g + ) as well as some Werner bands 
(C'llj^ <— X'E g + ), which allows for an accurate measurement 
of the H2 column densities in each component and in different 
rotational levels. A portion of the UVES spectrum covering the 
H2 Lyman (1-0) band is shown on Fig. [8] while the full veloc- 
ity plots for different rotational levels are shown on Figs. Q~6] to 
I2TI The measured column densities and corresponding excitation 
temperatures are given in Table |6j. In the following, we refer to 
these components as #1, #2 and #3 from the bluest to the reddest. 

Component #1 has similar width in all rotational levels and 
requires a single Doppler parameter to describe the profiles of 
the H2 J=0 to J=3 absorption lines. The Doppler parameter (b ~ 
3 km s~') is likely to be dominated by turbulent motions. 

Component #2 presents a broadening of the profiles with 
increasing rotational level and requires different Doppler pa- 
rameters. This beha viour has already been observed in the 
Galactic ISM (e.g. iJenkins & PeimberlL 11997b lLacour et all 
2005b) but also i n high redshift Damped Lyman-a' systems 
dNoterdaeme et all l2007al) . Doppler parameters can be mea- 
sured accurately even when significantly smaller than the spec- 
tral resolution thanks to the presence of numerous transitions 
with different oscillator strengths. However, the measurement 
of b in the first rotational level (J = 0) remains difficult due to 
the small number of unblended absorption lines (see Table [6). 
However, there are enough transitions from the J = and 1 lev- 
els together to ascertain the fact that the Doppler parameter of 
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Fig. 4. Portions of the J 1237+0647 UVES spectrum around Ci 
absorption lines. The short vertical marks in each panel represent 
the positions of absorption lines from the ground state ( 3 PoX 
the first excited level ( 3 Pi) and the second excited level CP2) 
from top to bottom, respectively. The component where CO is 
detected (at z a b s = 2.68957) is marked in red. 
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Fig. 5. Fit to S 1 absorption lines. Results of the best model fit 
using two components is overplotted (xl = 0.97). 
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is reached for logjV(S l) ) = 13.24, b = 0.91 kms" 1 . 
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Fig. 7. CI 1/11347 absorption at z abs = 2.68960. A single compo- 
nent fit is overplotted. 



the transitions from these two levels is very small, of the order 
of b ~ 1 km s _1 . This is consistent with thermal excitation with a 
temperature of ~ 120K, which is similar to what is measured 
from 7oi . 

The reddest component (#3) is very strong, with damping 
wings seen for rotational levels up to J = 3, allowing for accu- 
rate measurement of the column densities. Non-saturated lines 
for J = 4 and 5 also allow for accurate measurements in these 
rotational levels. Component #3 is particularly interesting not 
only because of its large H2 column density, but because it also 
contains deuterated molecular hydrogen, carbon monoxide as 
well as neutral sulphur and neutral chlorine, all of which have 
been very rarely detected at high redshift. Since the column den- 
sity of H2 is large, the J = and J = 1 levels are self-shielded 
and the collisional time-scale is much shorter than the photo- 
dissociation time-scale. Therefore, the Af(J=l)/Af(J=0) ratio is 
maintained at the Boltzmann equilibrium value. This means that 
the measurement of the kinetic temperature from 7oi is robust. 
We measure To\ = 7ki n ~ 1 10 K, which is similar to the temper- 
ature in the local interstellar medium (T^ ~ 80 K: ISavage et al.L 
fl977t) . 

We measure a total column density logA^(H2) = 19.21^*J-|| 
in the sub-DLA system with far the most important contribu- 
tion coming from the CO-bearing component (component #3). 
This corresponds to a large mean molecular fraction (fm) = 
2A^(H 2 )/(2Af(H 2 ) + N{H )) = 0.24+°^. However, the centre of 
the H 1 Lyman-a absorption line is clearly shifted from compo- 
nent #3 by more than 50 km s _1 (see Fig. [2j. This means that the 
amount of atomic hydrogen in the CO-bearing cloud is much 
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Table 6. H2 column densities and excitation temperatures 



Table 7. HD column densities 



component 

z abs , v (km s~') 




logA'(H 2 ,J) fl 
(cm" 2 ) 


b 

(kms-») 


To-, 
(K) 


#1, z=2.68801, A v=- 
J=0 


127 


16 28 +010 
15.51±0.05 


3.3±0.2 




J=l 
3=2 




16.08±0.10 
15.39±0.20 


>, 


193 + if 

v -54 

271+ i20 

-63 


1=3 




15.13+0.05 


» 


26U)« 


#2, z=2.68868, Av= 


-73 


17 62 +0 0S 






J=0 
J=l 




16 94 + ° 2( ! 
u J -0.44 

17.51+0.05 


4 + °? 
-0.1 

1.2+0.1 


> 117 


1=2 




15.25±0.03 


3.3±0.2 


93!f 


1=3 




14.89±0.02 


4.9±0.2 


131^> 


#3, z=2.68955, A v= 
J=0 


=-2 


19 20 +0 13 
18.65±0.20 


6.0±0.1 




J=l 
1=2 
1=3 
J=4 




18.92±0.10 
18.18+0.10 
18.21+0.10 
15.43±0.10 


7.9±0.1 


108^ 
190+g 
252 +52 
177 +14 


J=5 




14.95±0.05 




213 +u 



Notes. <a) The first line for each component gives the total H9 column 
density in that component. 



smaller than log N(H°) = 20 and the value given above should 
be considered as a lower limit on the actual molecular fraction in 
the CO-bearing cloud (i.e. /h 2 > 1/4)- As noticed in the previ- 
ous section from the presence of Cl° absorption associated with 
#3, the molecular fraction in the component #3 is probably not 
far from unity. 

3.4. HD and the D/H ratio 

Several Lyman-band HD lines from the first two rotational levels 
are detected in the UVES spectrum (see Fig. |9). Unfortunately, 
J = 1 lines are either severely blended or redshifted in regions 
of bad signal-to-noise ratio and we can only derive an upper 
limit on the column density for this rotational level. We mea- 
sure the HD column density in the J=0 rotational levelx from 
fitting the HD L5R0 and HD L8R0 lines which are unblended. 
Measurements are summarised in Table [7] 

We measure JV(HD) /2JV(H 2 ) = 0.95 x 10~ 5 . This is about 
10 times higher than what is measured in the Galactic ISM 
dLacour et all l2005al) for / = 0.24. Sin ce this ratio is known 
to increase with the molecular fraction dLacour et a [2005a) 
HD/2H2 might provide a lower limit on D/H for fm < 1 . If, as 
discussed previously, the actual molecular fraction in the HD- 
bearing cloud towards J 1237+0647 is higher than 0.25, both 
HD and H 2 could be self-shielded and HD/2H 2 ~ D°/H°. The 
value we obtain is then consistent with the D°/H° ratio mea- 
sured in the Galactic disc dLinskv et alll2006h . This corresponds 
to an astration factor of ~3 when comparing to the primor- 
dial value as measured in low-metallicity absorption systems 
(D/H = 2.82 ±0.2xlO" 5 :lPettini et al]l2008[|lvanchik et alJ[2O10L 
see however ISrianand et alj|2010l) or derived from the baryon 
density parameter ( Steigman, 2003). 



component 


logA'(HDJ) (cm" 2 ) 


Z>(kms-') 


z abs =2.68956, Av=-1 kms" 1 






J=0 


14.48+0.05 


4.5+0.2 


J=l 


< 13.60 





Note that all five high redshift HD detections to date yield 
relatively large D/H values despite significant metal enrichment: 
N(HD)/2N(H 2 ) = 1.5 x 10~ 5 , 3.6 x 10 ~ 5 7.9 x 10~ s and 
1.6xl0 ~ 5 towards resp ectively, J 1439+1 1 17 jNoterdaeme et all 
I2008bh . Q 1 232+082 rilvanchik et all 1201 Oh . J 2123-0500 and 
FJ 0812+32 (fTumlinson et all 1201 Oh . Since deuterium is easily 
destroyed as interstellar gas is cycled through stars, large deu- 
terium abundances are difficult to reproduce with closed-box 
models. However, these are well explained by models includ- 
ing infall of primordial gas (e.g. iProdanovic & Fieldsl 2008). 
If the velocity-metallicity correlation found by Ledoux et all 
(2006a) is the consequence of an underlying mass-metallicity 
relation, then we can expect that a high astration of deuterium 
in high metallicity systems is roughly compensated by a strong 
infall of primordial ma terial onto massive galaxies. However, 
Tumlinson et al] d2010t) noted that HD/2H 2 ratios in high-z ab- 
sorption systems lie in a narrow range well above the value mea- 
sured in the Galaxy while these systems present a large diver- 
sity in terms of metallicities and molecular fractions. This puz- 
zling behaviour led them to conclude that it could be premature 
to use the HD/2H 2 ratio to derive Q/,, given our actual under- 
standing of interstellar chemistry. In addition, we note that in 
the case of QSO absorbers, we only have access to the proper- 
ties of the gas (metallicities, molecular fractions) averaged over 
the line of sight. These may not be representative of the ac- 
tual chemical abundances in the HD-bearing cloud. Indeed, only 
the total 2V(H°) can usually be measured and the metal compo- 
nents are blended into a smooth absorption profile. It is therefore 
necessary to be careful and to study each system in detail (e.g. 
Balashev et al., submitted) to derive the local chemical and phys- 
ical conditions in the cloud. 



3.5. Carbon monoxide 

Absorptions from eight A l U(V) <- X'S + (v=0) bands of CO 
(from AX (0-0) to AX (7-0)), the C'Z^O) <- X'i: + (v=0) 
Rydberg band and the d 3 A(v'-5) <— X l I. + (v-Q) inter-band sys- 
tem are detected at z- d bs = 2.68957 towards J 1237+0647, see 
Fig. [10] Each band produces a complex absorption profile which 
results from superposition of absorption lines from different ro- 
tational levels in the P and R branches. The resolving power 
(R ~ 50000) and the signal-to-noise ratio (SNR ~ 28) of the 
UVES spectrum are high enough to individually measure the 
column densities in rotational levels up to J = 3. In addition, the 
J = 4 level is probably detected in the P branch of CX (0-0). We 
use the AX and dX bands that fall outside of the Lyman-a forest 
to measure the column densities in rotational levels up to J = 3. 
The AX (3-0) band is affected by a spike likely due to a cos- 
mic ray impact at the position of the R branch and this region is 
therefore not considered when fitting the profile. In addition, we 
use the CX(Q-Q) band up to J — 4. This band is the strongest one 
available but the corresponding rest wavelength (1088 A) makes 
it redshifted in the Ly-a forest. Fortunately, only the R branch is 
blended whereas the P branch is free from any blend. Moreover, 
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Fig. 8. Portion of the UVES spectrum of J 1237+0647 covering the (1-0) Lyman band of H2. The labels indicate the branches ('R', 
'P' for AJ = — 1, +1, respectively) and the rotational levels of the lower states. Absorptions from different components are indicated 
using different label colours (#1: blue, #2: green, #3: red). 
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Fig. 9. Fit to HD J=0 lines. The blue profile is the contribution 
from HD while the red profile represents H2. The velocity of the 
H2 detected components are indicated by vertical dotted lines. 



the CX rotational levels are well separated at the UVES spectral 
resolution. 

The results of the fit are presented in Table [8] Two errors 
are quoted in this table for the column densities: the first one 
is the rms error on the fit, while the second one reflects the un- 
certainties resulting from the continuum placement. The latter 
uncertainties were estimated by changing the normalisation by 
plus or minus 0.5 <x (i.e. about + 2% for SNR = 28) around the 
best continuum fit. The total CO column density we derive is 



logN(CO) = 14.17+0.09, implying a high CO to H 2 ratio of 
N(CO)/N(U 2 ) = IP" 5 - This is typic al of what is seen in translu- 
cent clouds (see Burgh e t al.Ll2010h . 



Table 8. CO column densities and excitation temperatures 



component 


logiV(CO,J) 


b (km s ') 


To-," (K) 


z abs = 2.68957 


14.17±0.09 c 


0.9+0.1 




J=0 


13.53±0.04±0.08 






J=l 


13.77±0.02±0.07 




10 i +4 - 3 

-2.0 


J=2 


13.54±0.03±0.06 




10 - 5 -n 


J=3 


13.21±0.04±0.09 




12 4 + ' ° 


J=4 


12.64+0.16+0.11 




13 +1 9 



Notes. ia) Quoted errors on column densities are respectively errors 
from fitting the Voigt profiles and errors due to continuum placement. 
The latter were estimated by varying the continuum by ±0.5<r. ib) The 
errors on T _j represent the extremum values for the different sets of 
continuum. (c) Total CO column density. 



In Fig. Q~T| we show the excitation diagram of CO. It is 
clear that the population of the first three rotational levels can 
be reproduced with a single excitation temperature. We mea- 
sure this excitation temperature by performing a linear fit of 
log N(CO,i)/gs vs the energy of the levels (Eqj)- The fit and 
1 cr range on Fig. QT] corresponds to the best fit continuum. In 
order to estimate the effect of the continuum placement, we re- 
peat the linear fit for each set of continua and take the extrema 
as representative of the range of possible values for r ex (CO). 
This gives T ex — 10.5^ g. Note that the effect of the continuum 
placement is mainly a change in the total CO column density, 
while little change on the slope of the linear fit (i.e. the exci- 
tation temperature). The CO excitation temperature is well be- 
low the kinetic temperature of the gas. This means that the gas 
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'Q') and the rotational level of the transition. 



density is small enough so that radiative processes are likely to 
dominate the downwards cascad e, as predicted in d iffuse molec- 
ular and translucent clouds (e.g. lWarin et all fl996h . Indeed, the 
population ratios of the neutral carbon fine-structure levels, in- 
dicate a volumic density of the order of n H o ~ 50 crrT 3 , well 
below the critical density at which the collisional de-excitation 
rate of CO(J=1-0) equals that of the sp ontaneous emission 
(« crit ~ 1000 cirT 3 ; ISnow & McCall 12006). Indeed, in terms of 
density and molecular fraction, the CO-b earing system p r esented 
here is very similar to that presented in ISrianand et alj (|2008b) 
where we concluded that collisional excitation of CO is negligi- 
ble. It is important to note however, that the fine-structure levels 
of C° only give the average volumic density. The actual local vo- 
lumic density in the CO-bearing cloud could be higher. A small 
shift (~ 1 km s~') is measured between the strongest C i feature 
and the CO component. This may indicate that the two species 
are not completely co-spatial. 

From the radiative code RADEX (Ivan der Tak et al ., 2007), 
we expect the excitation temperature of CO to be about one 
degree larger than the expected temperature of the Cosmic 
Microwave Background (CMB) radiation (7cmb(z=2.69) = 
10.05 K) as soon as the collision partner (H°, H2 and He) density 
is larger than 50 cm 4 . This explains that the excitation temper- 
ature we measure is slightly higher than what is expected from 
excitation by the CMB radiation alone. 




E 0J (K) 

Fig. 11. Excitation diagram of CO rotational levels. Errors on 
the column densities from fitting the lines are represented by the 
small red error bars while the long black error bars take into ac- 
count the uncertainty in the continuum placement. The plain line 
represents the linear regression fit using J = to J = 2 measure- 
ments and the best continuum fit. 1 <x errors are represented by 
dashed lines. 



In Fig. Q2] we compare the excitation temperature of CO at 
high redshift with that in the local Universe. In the local ISM, 
the temperature is seen a few degrees above Tqmb at low CO 
column densities and rises for column densities above N(CO)= 
2 x 10 15 cm 2 . This is due to the incre ased importance of pho- 
ton trapping at larger column densities (Burgh et al., 2007). The 
values observed at high redshift are significantly higher than the 
local ones, despite similar N(CO) and kinetic temperatures. This 
clearly means that the main physical difference between high 
redshift and local lines of sight is the higher CMB temperature 
at high redshift. This provides a strong positive test to the hot 
Big-Bang theory. Another consequence of Fig. [12] is that only 
CO-bearing systems with logTY(CO) < 15 - for which there is 
no correlation between N(CO) and r ex (CO) - are good places 
where to measure the evolution of 7cmb with cosmic time. 

Interestingly, although the differences are small and within 
errors, we measure a systematic trend, T04 ^ T03 ^ 7o2 55 Toi, 
regardless of the exact continuum placement. This indicates that 
while CMB photons dominate the rotational excitation of CO, 
other mechanisms are at play. We fail to reproduce the increas- 
ing temperature with increasing rotational level with RADEX. 
However, such behaviour has already been notice d in th e lo- 
cal ISM (ISonnentruckeret all l2007t ISheffer et all 120081) and 
could be explained by the select ive self-shielding of low ro- 
tational lines for logN(CO)>14 dWarin et all 1 19961) . The self- 
shielding of far-UV Rydberg bands of CO (those relevant to the 
photo-destructi on process) could b e more effective than previ- 
ously thought ( She ffer et all 120031) . In addition, the presence of 
H2 lines in the same spectral region can contribute to an effec- 
tive shielding of CO lines. Finally, radiative pumping from CO 
emission lines due to nearby dense molecular clouds could con- 
tribute to_gogulate_the higher rotational levels in the absorbing 
cloud dWannier et all 1 19971) . If the increasing temperature with 
increasing rotational level is physical, then Tqi = lO.l^n ^ 
could represent better the excitation by the CMB alone. 
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Fig. 12. Excitation temperature of CO as a function of the to- 
tal CO col umn density. Black error bars are measurements at 
z = from Burgh et al. (2007). The red filled circle with error 
bars is our measurement at z — 2.7 towards J 1237+0647 while 
the blue squar e represents the mea surement at z = 2.4 towards 
J 1439+ 1117 (ISrianand et all l2008bl) . Only a large range on T ex 
(6-16 K) could be determined for the system at z = 1.64 towards 
SDSS J160457.50+220300.5 dNoterdaeme et all l2009ah . and is 
therefore not represented in the figure. The vertical dotted line is 
indicative of a change of behaviour: below this limit there is no 
correlation between N(CO) and r ex (CO). 

4. The nature of the absorbing cloud 

4.1. Summary of the physical properties in the CO 
component 

In the previous sections, we have derived physical proper- 
ties of the gas associated to the molecular absorptions seen at 
z abs ~ 2.68957 towards J 1237+0647. 

From the analysis of Cl°, we derived that the molecular frac- 
tion, 2jV(H 2 )/(2jV(H 2 )+jV(H )), in the CO component is larger 
than 1/4 for a super-solar metallicity: Z(Zn,S) = +0.34,+0.15. 
From the populations of the C ground-state fine structure levels, 
we found that the particle density is of the order of ~50 cm -3 . 
The analysis of the H 2 rotational levels yields a kinetic tempera- 
ture of ~100 K and CO is mainly excited by radiative interaction 
with the CMBR. 

We can have an indication of the electronic density in the 
cloud thanks to the S°/S + ratio. Assuming that the mean ratio 
can be derived using the column densities in the strongest com- 
ponents we measure: log N(S°)/N(S + ) = -1.72. The electronic 
density, « e , is derived from the ionisation equilibrium between 
the two species 



r«(s°) = ff«e«(s + ), 



(i) 



where T is the photoionisation rate of S° and a the combina- 
tion rate of S + . Taking the ratio in diffuse ga s of the Galaxy 



(T la ~ 95 cm 3 , Pequignot & Aldrovandilfl9 86: see also table 8 



of lWeltv et aTI l999b) we derive n, 



1.85 (r/ff/95) cm" 3 . Note 
that this is an average value in the strongest metal component. 
Since S and S + are not co-spatial -as indicated by the different 
Doppler-parameters- the electron density derived here should be 
considered as an upper limit in the molecular gas. 




-0.2 0.0 0.2 0.4 



Fig. 13. Distribution of SDSS photometric g - r values for a sam- 
ple of 650 non-BAL QSOs at 2.681 < z < 2.881. The arrow 
indicates the position of J 1237+0647. 



We could apply the same reasoning to carbon but we cannot 
measure the C + column density as C H/11334 is highly saturated. 
In turn we can derive N(C + ) from the sulphur electronic density 
measurement. We find 



jy(c + ) _ r(c°)/g(c + ) jv(c°) 

N(S + ) ~ F(S°)/a(S + ) N(S°) 

Using F(C°)/a(C + ) = 32 cirr 3 yields logN(C + )/N(S + ) 
which is close to the solar ratio. 



(2) 



1.26 



4.2. Dust content 

Fig-EJshows the distribution of g - r colours for 650 non-BAL 
quasars with redshifts similar to J 1237+0647 (Az = 0.1). The 
median g — r value of the distribution is 0. 16 with a standard de- 
viation of 0. 1 1 . This dispersion reflects the color variation from 
one QSO to the other and is not due uncertainties in the SDSS 
photometric measurements (which are better than 0.03 mag). 
This implies that the colour excess of J 1237+0647 (g—r = 0.59) 
is significant at the 3.9 cr level. This significance increases to 
4.9 cr if we consider a Gaussian fit to the distribution. Indeed, 
there is a tail towards large colour excesses which shows the 
existence of a population of reddened quasars. J 1237+0647 is 
among the reddest 1.5% quasars having g - r > 0.5. 

The observed flux-calibrated SDSS spectrum of 
J 123 7+0647 is matched by the QSO composite spectrum 
from lVanden Berk et all d2001l) reddened by a SMC-like extinc- 
tion law with E(B-V) = 0.05+0.01. The spectrum together with 
the fit is shown in Fig. [14] Note that since we are only interested 
in the slope of the continuum of J 1237+0647 and the strength 
of emission lines varies strongly from one quasar to the other, 
we do not include emission line regions in our fit. Absorption 
lines were also ignored during the fit as well as the wavelength 
range bluewards of the Ly-a emission. 

In order to further estimate the probability that the above 
reddening might be due to a peculia r intrinsic QSO shape , 
we used a technique described in (ISrianand et all |2008a; 
iNoterdaeme et alll2009al) : We repeated the spectral slope fitting 



11 



P. Noterdaeme et al.: A translucent interstellar cloud at z = 2.69 



assuming an absorber at z = 2.69 for a control sample of 82 non- 
BAL QSOs with similar emission redshifts (2.76 < z em < 2.80) 
and spectra with /-band S/N ratios larger than 5. We find that the 
continuum slope of J 1237+0647 deviates at the 98% confidence 
level from the mean slope of other quasars. In the following, we 
will therefore consider that the colour excess of J 1237+0647 is 
due to the presence of dust at z a b s = 2.69. 

Our best-fit model predicts J = 16.9, H = 16.3, K = 15.5 
(with typical errors of 0.2 mag), while the observed Two 
Micron All Sky Survey (2MASS) magnitudes are J = 16.8+0.2, 
H = 15.7+0.2, and K > 15.0. The agreement although not per- 
fect in the H-band is reasonable. Indeed, the 2MASS magni- 
tudes come from Point Source Reject Table for objects with 
very low SNR. These measurements are known to suffer from 
flux overestimation which can easily explain the discrepancy 
between predicted and measured H-band magnitudes. In addi- 
tion, the presence of the H-/3 emission line in the H-band in- 
creases the uncertainty of our estimate. This together with the 
fact that the 2MASS and SDSS observations were taken five 
years apart, makes our predicted magnitudes consistent with the 
2MASS data. Unfortunately, there are no SDSS measurements 
at different epochs for this object to monitor any variation in the 
QSO flux. 

The measured reddening, although significative, is 
marginally higher to what is seen in the ge neral popula- 
tion of DLAs (E(B-V) ~ 0.04. lEllisonetal.ll2005l) . Interestingly, 
the integrated extinction-to-gas ratio measured towards 
J 1237+0647, A(V)/7V(H°) = 1.5X10 -21 mag cm 2 is 20 times 
high er that the avera ge value for the SMC (7.5xl0~ 23 mag cm 2 ; 
iGordon et all l200l and about 50 times higher than the mean 
value measured in high redshift DLAs (2-4xl0 -23 mag cm 2 ; 
IVladilo et all 120081) . This means that would the H° column 
density have been higher, the extinction induced would have 
been so large that the QSO would have been missed by the 
SDSS target selection^- Note that the moderate extinction 
Ay = 0.14 in the rest-frame of the absorber, already produces an 
extinction of nearly 1 mag in the g-band. This supports further 
the possibility that current surveys can miss a large number of 
cold clouds dNoterdaeme et all 1200 9 a). If, as discussed before, 
a large fraction of the atomic hydrogen is not associated with 
the molecular component, then the extinction-to-gas ratio in 
the molecular component is even higher. The line-of-sight to 
J 1237+0647 probably passes through a relatively thin slab 
of dusty gas. This is also supported by the high depletion 
factors measured in the CO-bearing component ([Fe/Zn] — 2.0, 
[Ni/Zn]~-1.6, [Si/Zn]~-1.7, [Cr/Zn] < -1.6). Such abundance 
pattern is typical of what is seen in cold gas of the Galactic disk 
(IWeltv et all[T999ah . 

If we consider the extinction only, then the sightline studied 
in this paper is not translucent. Indeed, the historical definition 
of translucent correspon ds to 1 < Ay < 5. However, as noted 
bv lRachfordetal.1 £2002), even translucent sightlines can result 
from the concatenation of multiple diffuse clouds along the line- 
of-sight. This kind of scenario has already been advocated to 
explain the bimodal distribution in the logjV(H?) dis t ributio n 
for a given column density of dust dNoterdaeme et all [2008a). 
Therefore extinctio n may not be the bes t par ameter to define 
translu cent clouds. iRachford et al.l (120021) and ISnow & McCalll 
(2006) have proposed that the definition of translucent clouds 
should be based on the local properties of the gas rather than 



2 For logA'(H°) = 20.65 and same extinction-to-gas ratio, the pre- 
dicted magnitude already reaches the ( = 20.2 limit set by the SDSS-II 
collaboration for spectroscopy of high redshift quasars. 




4000 5000 6000 7000 _ 8000 9000 

Observed wavelength (A) 

Fig. 14. SDSS spectrum of J 1237+0647 (black). The orange 
dashe d spectrum is the unredden ed SDSS QSO composite spec- 
trum (IVanden Berk et all l200ll) . The red spectrum is the com- 
posite spectrum reddened by a SMC extinction law with E(B- 
V) = 0.053 (i.e. Ay = 0.14, in the rest-frame of the absorber). 
The fit was perfo rmed using an IDL code based on MPFIT 
dMarkwardlll2009l) . 



the integrat ed properties alon g the line-of-sight. Following this 
suggestion, Burgh et al. (2010) have presented a definition based 
on the abundance of hydrogen and carbon in molecular forms. 
These authors noted that with this definition, translucent clouds 
do not produce necessarily strong reddening. 

4.3. The molecular fractions of H and C 

In Fig. [15] we plot the ratios N(CO)/N(H 2 ) and N(CO)/N(C°) 
versus the hydrogen molecular fraction, fui, for the two systems 
towards J 1237+0647 (this work) and J 1439+1 1 17 (Srianand et 
al. 2008) together w ith measurements in the ISM of the Galaxy 
(Bur gh etal.L l2010). The latter authors define translucent clouds 
as clouds with CO)/H 2 > 10~ 6 and CO/C° > 1.0 for f m > 0.4. 
In the cloud toward J 1237+0647, we find N(CO)/N(H 2 ) = 10~ 5 
for (fm) - 0.24. However, as discussed previously, the hydrogen 
molecular fraction in the CO-bearing cloud is larger than this. In 
addition, the amount of carbon in CO molecules is about that in 
atomic formB We therefore conclude that the cloud in front of 
J 1237+0647 is indeed a translucent cloud, an ideal laboratory 
from probing astrochemistry at high redshift. 

5. Conclusion 

From our VLT survey for H 2 in DLAs, it appears that neutral 
carbon is generally observed in the same components tha t fea- 
ture H 2 dSrianand et all 120051; iNoterdaeme et all l2008al) . We 
therefore selected the rare SDSS lines of sight in which Ci ab- 
sorptions are present. From UVES follow-up observations, we 
have detected strong absorptions from H 2 , HD and CO along 
SDSSJ123714.60+064759.5. This is a beautiful and peculiar 
case where detailed analysis of the physical properties of the gas 
is possible. 

The H° column density is small, log N(B°) = 2 00+0.15 and 
corre sponds to what is usually called a sub-DLA dPeroux et all 
2002). Corresponding overall metallicity is super solar with 
[Zn/H] = +0.34 and [S/H] =+0.15 . The system features three 
H 2 components spanning ~125 km s _1 , the strongest of which, 

3 Note that we us e here log N(C°, 7 = 0, 1 , 2) = 14.39 ± 0.10 derived 
using f- val ues fromlJenkins & T ripp (2001) to enable comparison with 
the work bv lBurgh et alJ d2010h 
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Fig. 15. N(CO)/N(H 2 ) and N(CO)/N(C°) versus the molecular 
fraction, fm- The red filled circle with error bars is our mea- 
surement at z = 2.7 towards J 1237+0647 while the blue square 
represents the measure ment at z = 2.4 towards J 1439+1117 
(ISrianand et all 1200 8b). Other points represent measurements 
in the Galactic IS M with three outliers represented in grey 
(Bur gh et all 12010). These are peculiar systems with uncertain 
measurements by the same authors. The horizontal dotted line 
marks the limit between diffuse and translucent clouds, while 
the two vertical dotted lines mark the transition range between 
these two regimes. Note that the molecular fraction we indicates 
toward J 1237+0647 is a lower limit. In order to enable easy 
comparison wi th local values, we use N(C°) measured with the 
/-values from Jenkins & Trip3 d2001l) . 



with logAf(H2) = 19.20, does not coincide with the centre of 
the H i absorption. This means that the molecular fraction in this 
component is larger than the mean molecular fraction (fm) = 1/4 
in the system. 

From the populations of the low H2 rotational levels, we 
measure the kinetic temperature of the gas to be around 100 K in 
the strongest component, where HD and CO are also detected. 
The detection of S and C° implies that the gas is shielded from 
the surrounding far-UV radiation field. 



The relative populations of the C° fine structure levels yields 
an estimate of the average hydrogen density in the main com- 
ponent of about 50-60 cm 4 . At such densities, collisions are 
not frequent enough to dominate the rotational excitation of 
CO molecules and radiative processes are likely to determine 
the CO rotational populations. The excitation temperature we 
measure (T ex (CO)=10.5 K) is significantly larger than that mea- 
sured in the local ISM (r ex (CO) ~ 3.5 K) for similar CO col- 
umn densities, molecular fractions and kinetic temperatures. 
We show that the higher rotational excitation of CO towards 
J 1237+0647 results from the higher temperature of the Cosmic 
Microwave Background at the redshift of the absorber, 7cmb(z = 
2.69) = 10.05 K. This provides a strong positive test to the hot 
Big-Bang theory. 

Small velocity shifts (of the order of 1 km s _1 ) are observed 
between the different neutral and molecular species in the main 
component. For H2, the shift might appear more important but 
is rather due to the H2 component being a blend of several 
sub-components. Neutral chlorine i s likely a better indicator o f 
the strongest molecul ar component dSonnentrucker et ail l2006). 
ISrianand et all d2010h recently observed similar small velocity 
shifts between H2 and 21 -cm absorption in a high redshift DLA 
also testifying the presence of inhomogeneities of the ISM on 
very small scales. In turn, since the distribution of H° and O 
are expected to closely follow each other, the 1/11302 absorp- 
tion profile indicates that atomic hydrogen is distributed over the 
full ~400 km s _1 velocity range over which metals are observed. 
This explains the large velocity shift between the centroid of H 1 
and molecular lines. 

All this reinforces the view that the ISM is patchy, with 
small and dense molecular cloudlets (probably with an onion- 
like s tructure) immersed i n a warmer diffuse atomic medium (see 
e.g. |Petitjean et all 119921) . We can derive an upper limit on the 
size of the molecular-rich region along the line-of-sight by con- 
sidering that all H° is associated with the main velocity compo- 
nent of density «h = 50 cm 4 . The corresponding characteristic 
length is / = jV (H )/w H o - 0.6 pc . A lower- limit can be put us- 
ing E q. A7 of Uura & Yorkl (l978, see also ISonnentrucker et all 
2002) and the observed molecular and neutral chlorine fractions. 
This gives I > 0.05 pc. The molecular region of the system 
has therefore a very small size and hence small cross-section. 
It is therefore not surprising that detections of translucent clouds 
were elusive till now. Studying the frequency of CO absorbers 
would give an idea of the filling factor of the molecular-rich 
gas, but requires larger statistics. Small physical extents could 
yield partial covering of the background source by the cloud. 
This may happen in particular if some absorption lines are red- 
shifted on top of emission lines from the extend ed QSO broad 
line r egion, as seen in the case of Q 1232+082 (llvanchik et al.l 
120101: Balashev et al., submitted). 

Interestingly, the conclusion that the ISM at high redshift 
is made of small cold cloud lets immersed in w armer diffuse 
medium has been reached by iGupta et"aTI d2009t) while consid- 
ering the distribution of cold gas detected in 21 -cm absorption. 
Note that J 1237+0647 is detected in the radio domain by FIRST. 
However, the low radio flux (2.3 mJy) prevents any spectro- 
scopic study with current radio-telescopes. 

Several chemical reactions can take place in this cloud. 
Indeed, three molecular species are detected, while the pres- 
ence of neutral chlo rine suggests chemica l reactions involving 
HC1 and HC1 + (e.g. Dalga rno et all 11974 iNeufeld & WolfireL 
2009). The DLA system toward J 1237+0647 is therefore an ex- 
cellent candidate to target with future extremely large telescopes 
(ELTs) to detect other molecular species like C2, CH, OH and 
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study astrochemistry in the interstellar medium of high redshift 
galaxies. Given the expected attenuation of the quasar by CO- 
bearing clouds, X-shooter, with its high throughput and medium 
resolution, is the best instrument to survey carefully selected 
DLA/sub-DLA candidates for CO absorptions and build a sam- 
ple of molecular-rich clouds at high redshift, which may then be 
studied in details in the optical and sub-millimeter ranges with 
future facilities like ELTs and ALMA. 
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Fig. 16. Fit to H2(J=0) lines. The blue profile is the contribution 
from HD. 
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Fig.l7.FittoH 2 (J=l)lines. 
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